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Abstract 



We present a derivation of the probability of bremsstrahhmg from high- 
energy electrons in a screened Coulomb field using the quasiclassical op- 
r S^ . erator method. It is shown that recent Zakharov's criticism is completely 

groundless and comes from misunderstanding of the method. We confirm 
all our published results. 



1. Comment, recently published by Zakharov [|I|, shows that it is reasonable 
to discuss once more the basis of the quasiclassical operator (QO) method devel- 
oped by authors in |2|, |3|]. The method was presented in detail in the book || 
and later in the books and |J (the latter is the enlarged English translation 
of the book ||). Very clear description of the QO method is given in the popular 
textbook fTJ. 

This paper is divided into two parts. In the first part we discuss in frame 
of the QO method the derivation of the basic formulae describing the radiation 
process paying special attention to the important details of the approach. In 
the second part we analyze Zakharov's claims step by step and show that all 
his pretensions are erroneous because of complete misunderstanding of the QO 
method. 

In this paper we consider the applicability of the QO method to the problem 
of radiation from a ultrarelativistic particle at a potential scattering. In this 
process the emitted photon and the final electron are moving at a small angle to 
the initial electron momentum and the large angular momenta I ^> 1 contribute. 
In this situation the quasiclassical scattering theory is applicable. The very 
important result of the method is that recoil at radiation is incorporated into 
the theory in the universal form for any external field. There are two essentially 
different cases in application of the QO method. In the first case scattering in 
a potential can be considered in classical terms: phase shifts are large, there is 
a correspondence between the impact parameter and the momentum transfer. 
Therefore it is possible to use the version of the QO method where one can 
instead of operators substitute classical variables in the coordinate space in an 
expression for probability of the process (see, @, Sees. 9-15 in ||], Sees. 2-5 in ||, 
HI, Sec. 90 in 0). In the second case, the process of scattering can't be described 
in classical terms. However, for scattering of ultrarelativistic particles where large 
angular momentum contributes, one can use the quasiclassical approximation for 
description of scattering including the situation where the phase shifts are small 
(see §, Secs.9, 16-20 in [g, Secs.2,7 in §, @, Sec.96 in §). Precisely this 
formulation of the method must be applied for consideration of radiation from 
ultrarelativistic particles at scattering on atoms in a media. In our recent papers 
devoted to the theory of Landau-Pomeranchuk-Migdal (LPM) effect [|J-|jrTJ we 
use the mentioned formulation of the QO method. Therefore we will analyze it 
in detail below. 

2. Because the LPM effect can be observed at very high energies, we can 
consider the case of complete screening, so that (we employ units such that 
h = c = l) 

a s < = — it 1 = k, (1) 

q min ujm 2 

where a s is the screening radius (a s ~ 111Z _1 / 3 A C , A c = l/m), Z is the charge of 
a nucleus, q m in is the minimal momentum transfer which is longitudinal (with 
respect to the momentum of the initial electron p), e(u) is the energy of the 



initial electron (the emitted photon), m is the electron mass, Iq is the radiation 
formation length for a small angle scattering on an isolated atom. Note that 
in frame of the QO method the radiation problem is solved for the case of an 
arbitrary screening , see Sec. 18 in ||. The impact parameters g, contributing 
into the scattering cross section, are small comparing the formation length (g < 
a s ^ Iq) i n a screened Coulomb potential. This means that the scattering of 
the ultrarelativistic particles (the virtual electron is close to the mass shell) 
takes place independently of radiation process (see, Sec. 7 in [|], ||, Sec. 96 in 
[0]). Thus, we can present the cross section of radiation as a product of the 
probability of emission of a photon with the momentum k at given momentum 
transfer q^ (q±p = 0), and the cross section of scattering da(q±) of a particle 
with the same momentum transfer qjj 

do 1 = W 7 (qj_, k)d 3 kda(q 1 _). (2) 

We show below that in frame of the QO method the probability of radiation 
W 7 (qj_,k) is given by the classical trajectory of a particle in "the form of an 
angle" in the momentum space 

p(t)=#(-*)P + 0(*)(P + <U)> (3) 

while the cross section do~(c{±) should be taken in the eikonal form. 

3. The probability of photon emission from an electron in frame of the QO 
method has the form (see, Eqs.(16.7)-(16.12) in g, Eqs.(7.1)-(7.12) in [|, §, 
Eqs.(96.1)-(96.8) in 0) 

dw =< i\M + M\i >, (4) 

where 

M = — %= [ R(t) exp i f ^-dt' 



o 



(5) 



here R(t) = R(p(t)), -R(p) is the matrix element for the free particles depending 
on the electron spin, kp = ooTi — kp, 7i = y/p' 2 + m' 2 , \i > is the state vector of 
the initial particle at the time t = 0, and p(t) is the operator of momentum in 
the Heisenberg picture: 

p(t) = exp(-iHt)pexp(iHt), H = H + V(r), (6) 

where V(r) is the potential of an atom. 
We present the evolution operator as 

exp(-iHt) = exp(-iHt)N(t), N(t) = exp(iHt) exp(-i(ft + V)t). (7) 

Differentiating the last expression over the time we obtain 

^-t- = -i exp{im)V{r) exp(-i(W + V)t) = -iV{r + vt)N(t),v = ^-. (8) 
dt rl 



The solution of this differential equation for the initial condition N(0) = 1 is 



N(t) = Texp 



-i f V(r + vt')dt' 



(9) 



where T is the operator of the chronological product. Bearing in mind that the 

commutator 

i 

[ r i> v 3\ = ft( S H- v i v j)> ( 10 ) 

one can drop operator T of the chronological product within relativistic accuracy 
(i.e. with accuracy up to terms ~ I/7) and present operator N(t) as 



N(t) — exp 



T 

: fv(g,z + t')dt' 



(11) 



where the axis z is directed along the momentum of the initial particle. If the 
formation time of radiation is much longer than the characteristic time of the 
scattering, one can present the dependence of the operator p(£) on the time in 
Eq.QD as 

p(t) = tf(-t)p(-oo) + 0(t)p(oo), p(±oo) = iV + (±cx))piV(±oo). (12) 

It should be pointed out that in the case when the scattering process is of non- 
classical character the operators p(— 00) and p(oo) are noncommutative among 
themselves and, generally speaking, one can't neglect their commutator. Let us 
mention also that using Eq.flTTD one find 



z oo 

P_l(-oo) ~ p± + J V b V{q, z)dz\ p ± (oo) - P± - / V g V(g, z')dz' . (13) 

—00 2 

These approximate expressions are written in the classical form when the entering 
operators commutate among themselves. Substituting the "trajectory" ( |T2"D in 
Eq.flD we obtain (this is Eq.(7.7) in [|, §) 



M 



ie{£ — us) 



i?(p(oo)) i?(p(-oo)) 



kp(oo) 



kp(— 00) 



(14) 



The state \i > in Eq.(^j) denotes the wave function in the configuration space 
in the time t — 0. However, in the scattering problem under consideration 
(motion of a particle in a local potential), the wave function are defined at 
t — > ±00. Let us consider the interconnection between these wave functions. By 
definition of the state evolution, we have 



lim exp(— iHi)\ 00 >= lim exp(— iHt)\0 >, 



|0 >= iV + (oo)|oo >= I out >~ exp 



V(g,z')dz' 



CXD > 



(15) 



For the state evolving from t — > — oo we have 



|0 >= N + (-oo) 



oo >= in >~ exp 



V(Q,z')dz' 



oo > . 



(16) 



As a rule, the state at t — > ±00 are the plane waves, i.e. they are the eigenstates 
of the momentum operator p. However, one can use also another definition 
of these states, e.g. wave packets. Such definition was used in description of 
the bremsstrahlung at the collision of two high-energy beams with the restricted 
transverse dimensions [ 12| . If there is a summation over the set of states, i.e. they 
are the intermediate states, choice of the states is determined by a convenience of 



calculation. This circumstance was used in 16 for the derivation of connection 



between spectra of the bremsstrahlung and the pair creation in the theory exact 
in Za. 

4. It follows from the above analysis that for the derivation of the differential 
cross section of the bremsstrahlung in Eqs.(|2]), ([|), it is necessary to insert the 
projection operator \f >< f\, which corresponds "out" states of the plane wave 
with the momentum pj, between the operators M + and M, and to take "in" 
states with the momentum p^ as the states \i >. So, the initial state is the 
eigenvector of the operator p(— 00) and the final state is the eigenvector of the 
operator p(oo): 



p(-oo)|i >= iV + (-oo)piV(-oo)iV + (-oo)|pi >= p { \i >, 
p(oo)|/ >= N + (oo)pN(oo)N+(oo)\p f >= Pf \f > . 



Using (J17D we deduce for the matrix element of the operator M (|T 

~R(p f ) R(p 



M 



ie\e — lu) 
2n\/uj 



kp f kpi 



<f\i> 



(17) 



(18) 



Using the definitions of "in" and "out" states Eqs.(|l5|),(|T6|) we have 



< f\i >=< P/ |AT(oo)iV + (-oo)|p i >=< p/|Texp 



-i f V(r + vt)dt 



Pi > 



t< p/| exp 



00 00 



-1 J V(r + vt)dt ~\ J dt 2 J dtx\y{h),V{t x )\ 



-00 — 00 



\Pi >, (19) 



where V(ti) means V(r + vti), etc. We retain here only the first two terms in 
the decomposition of the T-product being restricted to the commutator 



[V(t 2 ), V{t x )} = --(t 2 - t 1 )V ± V(h)V ± V(t 1 ) 



(20) 



The first term in the decomposition of the T-product in Eq. (|T^) gives the scatter- 
ing amplitude in the eikonal approximation, the second term gives the correction 
to the eikonal approximation: 

< f\i >= J d 2 gexp[iq ± Q + ix(Q)]27i5(p f i\ -pq), (21) 



where 



X(Q) = Xo(ff) + Xi(q), Xo(q) = ~ I V(q, z)dz 

— oo 
oo oo 

Xi{Q) = Ye I dZ2 ! ^iV^^Vx^x)^ - z x )d{z 2 - z x ). (22) 

— oo — oo 

In the centrally symmetric potential we have 

g 2 J dV 2 {^Q r Tz I ) 



xo(e) = - j v(y/? + z*)dz, X i(e) = -j j % dz ( 23 ) 



-oo 



In the Coulomb field xi = {Za) 2 7r/(2ge). The main contribution into corrections 
to the probability of the bremsstrahlung connected with the phase xi gives the 
impact parameters g ~ A c = 1/m, so that xi ~ (Z a ) 2 /l- The corresponding 
corrections ~ I/7 to the bremsstrahlung cross section were calculated by authors 
in 
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We introduce now the notations p« = p, pf = p' + k where p' is the momen- 
tum of electron after photon emission. Then neglecting the terms of the order 
of q\\ in the argument of the 5- function in Eq. ( pT| ) we have in the region q± 3> q\\ 
which contributes for the potential considered 



5{p n - Pl \\)~5{e' + u-e), e' = ^/p' 2 + m 2 . (24) 

Using this relation one can express the propagator kpf through the propagator 
kp' (up to terms of the order I/7 2 ) 



kp f = wy (p 7 + k) 2 + m 2 - k(p' + k) 

l / kp'\ e' 

= u)\Je 2 — 2kp' + kp' — ue ~ uie 1 1 — + kp' — ue = —kp'. (25) 



Our final result for the differential probability of radiation in Eq. (|2|) is therefore 

(26) 



a 1 

iy 7 (q ± ,k) = — — - 

{2ir) z to 



sR(p' + k) e'R(p) 
kp' kp 



The explicit expression for the probability W 7 (q_i_, k) with regard for polarization 
and spin effects are given in |4[. 

5. In the above analysis we traced in detail the transition to the expressions 
calculated on the trajectory of particle in "the form of an angle" in the momen- 
tum space. Precisely these definite trajectories (taking into account the recoil at 
the radiation) determine the probability of photon emission 



dw 



a 



d 3 k 



(2tt) 2 lu 

3 i. - P 



dt 2 I dtxR*(t 2 )R(ti)ex.p 



t-2 



re 



kv(t)dt 



W 7 d 6 k, v 



;i,v), 



e — uj. 



(27) 



When the projectile is moving in a medium it scatters on atoms. In this case 
the probability should be averaged over all possible trajectories with the weight 
function defined by the cross section dcr(q_i_). In the presence of a macroscopic 
external field there is the systematic variation of the mean velocity of a projectile. 
The probability of radiation in this case is given by the same Eq . (|27|) . Both, 
the systematic variation of the mean velocity of the projectile as well as the 
fluctuations of it velocity due to scattering with respect to the mean value are 
taken into account with the aid of the distribution function. Thus, when the 
projectile is moving in a medium the probability of photon emission per unit 
time has the form (see, Eq.(20.31) in @] and Eq.(2.4) in fTJ) 



dw 



' dw\ a d 3 k 

~dt] ~ (2vr) 2 u 



oo 

Re / dr exp I —i—tur J / d 3 v / d 



V 



f d 3 r f d^'Ci'd, tf^Fiir, v, t)F f (r', v', r; r, v) exp 



i— k(r' — r) 



(28) 



We consider here the case of an infinitely thick target. The distribution functions 
Fi and F f satisfy the kinetic equation (see, Eq.(20.32) in j| and Eq.(2.5) in 0)): 

dF(r,v,t) dF(r,v,t) dF(r,v,t) 



n 



dt dr dv 

Ja(vy)iF(ry,t)-F(v,v,t)}d 3 v', 



(29) 



where n is the number density of atoms in the medium, and <r(v, v') is the 
scattering cross section which in our case should be calculated in the eikonal 
approximation in the screened Coulomb potential. The combinations entering 
Eq.(|7|) are 



t-2 



tl 



kv(t)dt 



UJ 



kv(t))dt 



LOT 



k(r(i 2 ) - r(iO) 



LOT 



k(r' 



ti 



R*(t 2 )R(t 



1 



r. | 777, 



2e' 2 



j2m2 + (e> + e")M 



e 2 



\C(*,&). (30) 



Just this approach was used in papers |]S|-[IT| in the derivation of the principal 
equations in the theory of the LPM effect. 



Let us analyze specific Zakharov's claims in p. 4 of |l.| (given in quotation 
marks) . 

1. ".. Eq.(l) (this is our Eq . (p7|) ) derived neglecting the variation of the field 
acting on electron...". This statement is wrong. The Eq. (|27|) , is valid for 
any variation of the external field. Although we consider the case of slowly 
varying external field in the first derivation || , nevertheless for 7>1 (this 
is the only case we consider), this result can used for the Coulomb field (see 
Eqs.(10) and (11) in ||). Later we show that Eq. (|27|) , is valid for any field 
(see ||, 0) including the field of the plane wave where the dependence on 
both coordinate and time is essential. 

2. "..In the quantum regime the above correspondence between the scattering 
angle and the electron impact parameter in interaction with a medium 
consistent is lost. In this case one must take into account accurately the 
variation of the field acting on the electron in evaluating the radiation 
rate. Eq.(l) (obtained for a smooth field) does not make any sense in the 
quantum regime. It is well known that for a screened Coulomb potential 
the quasiclassical situation takes place at Za ^> 1. In real media, when 
Za < 1 we have essentially quantum regime. For this reason, it is evident 
that Eqs.(l) and (2) (Eq.(2) we will discuss below) are not justified for 
real media. The use of Eqs.(l) and (2) which are not valid in the quantum 
regime is the main conceptual error in the treatment of the LPM effect 
given in Refs.fll, 12]. (here this is Refs.§ and 0)." 

A few errors accumulate here. First of all it is evident that Zakharov is 



not familiar with classical papers of Bethe and Maximon [1J| (see also |fl£ 
flr§| ) where cross sections of the bremsstrahlung and the pair creation in 
the Coulomb field for high energy region are calculated with the use of the 
quasiclassical wave functions. The results obtained in the mentioned papers 
are valid for any value of Za, for low Z they coincide with the Bethe- Heitler 
cross section, for Za ~ 1 the cross sections contain additional term which 
is the Coulomb correction. We want to stress that our approach has the 
same region of applicability as the Bethe-Maximon theory. This theory is 
developed in frame of Schrodinger picture, while we are developed approach 
in frame of Heisenberg (operator) picture. Because of this, starting with 
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some stage of calculation our results coincide with those of Bethe and 



Maximon or Olsen, Maximon, and Wergeland |Ty[ . The case of a screened 
Coulomb field is considered in the last paper. 

Second, as it is shown above, our approach is basically quantum (quasiclas- 
sical) and only the final result can be in many cases, expressed in classical 
terms as it was explained above (and in @-@). Moreover, in the papers 
H, the case of the Coulomb field was considered. For calculation of the 
probability of radiation in this field for q ~ q m in one can't use the trajec- 
tory "in the form of an angle", the impact parameters g ~ Iq contribute 
(see Eq.([TJ)), and the cross section of radiation can't be presented in the 
form Eq.(0). But even in this case the Eqs.(f|) and (|j) are valid. The 
underlined claim is erroneous. In fact, the quantum- mechanical properties 
of scattering are relevant part of our theory. 

3. point 2 "The authors of Refs. [11,12] (this is Refs.[Q, |L4| here) compensate 



the incorrectness of Eq.(l) (this is Eq.(|27|) here) by another error in the 
procedure of averaging over the electron trajectories. According to their 
logic the QO expression (1) (this is Eq.(|27|) here) must be averaged over 
all possible classical trajectories (and author say this)." The last state- 
ment has nothing common with our method. As it was shown above the 
Eq. (|27| ) is valid in all cases when the motion of particle can be described by 
quasiclassical quantum mechanics. As was explained above, the averaging 
procedure appears only in the case when we have to consider the scatter- 
ing of the projectile in a medium. In this situation the radiation process is 
described by Eq. (|28|) . In a very definite form this statement is contained 
in [Q (the paragraph before Eq.(2.4)) 



Let us make the final outlook on the discussion. 

• The method developed in [0-@ and presented in detail in @]-[0 is basically 
quantum theory in Heisenberg picture in quasiclassical approximation. So, 
in use of the theory one have to proceed with operator calculus. In many 
case, but not always, the final result can be expressed in classical terms 
(also in a very definite form). Of course, the possibility to use classical form 
is a great advantage of the method since it simplifies essentially a solution 
of a particular problem and permits to solve problems which could not be 
solved by other methods. 



Zakharov [[I]] expressed pretensions to Eq.(2.1) of [13]. This formula is given 
as known with the reference to Eq.(9.27) of j|. So, to use this equation 
one have to read, at least, the mentioned Sec. 9 of this book. The problem 
was considered in detail also in books |5|], || which contain derivation 
of the mentioned formula. However these books were not mentioned in 
1 1 1 as well as known textbook KJ. Instead of this Zakharov invents his 
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own prescriptions which have nothing common with our method and then 
criticizes these prescriptions. See, e.g. p. 4, point 2 "...the distribution 
function... should satisfy the kinetic equation in which the collisional term 
contains the classical scattering cross section" . 

It is evident that there is essential difference between the classical theory 
and the quantum theory in quasiclassical approximation. Unfortunately, 
from Comment [ij it looks that Zakharov don't distinguish between these 
two. 

It is curious that any references on incorrect expressions in our publications 
are absent in Comment [I]]. 



It is rather strange to read that in 1996 Zakharov in [19| quite successfully 



derived formulae, in absence of an external field, contained in our paper 



14| and book HL Moreover, he explains what kind of cross section should 



be substituted (just the formulae used in the mentioned publications). 

It should be noted that although Zakharov reproduced some of our results 
many years later than they were obtained in original papers, he never refer 
to our publications before and only in Comment [[I]] there is citing of our 
publications. 



It's amusing that in list of criticized papers there is paper EH] devoted to 



quantum theory of the transition radiation where there is no scattering at 
all (the criticism devoted to description of scattering) and this once more 
shows the level of the pretensions. 
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